Influence of Phonon dimensionality on Electron Energy Relaxation 
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We studied experimentally the role of phonon dimensionality on electron-phonon (e-p) interaction 
in thin copper wires evaporated either on suspended silicon nitride membranes or on bulk substrates, 
at sub-Kelvin temperatures. The power emitted from electrons to phonons was measured using 
sensitive normal metal-insulator-superconductor (NIS) tunnel junction thermometers. Membrane 
thicknesses ranging from 30 nm to 750 nm were used to clearly see the onset of the effects of two- 
dimensional (2D) phonon system. We observed for the first time that a 2D phonon spectrum clearly 
changes the temperature dependence and strength of the e-p scattering rate, with the interaction 
becoming stronger at the lowest temperatures below ~ 0.5 K for the 30 nm membranes. 

PACS numbers: 63.22.+m, 63.20.Kr, 85.85. +j 



It is an established fact that at sub-Kelvin tempera- 
tures the thermal coupling between conduction electrons 
and the lattice becomes very weak This has signifi- 
cant implications for the operation of low-temperature 
detectors and coolers [2j, or for any solid-state sys- 
tems where dissipation and cooling are relevant. Low- 
temperature electron-phonon (e-p) interaction has been 
studied widely during the past decades, but mostly only 
for the case in which the phonons are fully three di- 
mensional (3D) d H H Q. However, due to signifi- 
cant advances in fabrication of thin suspended structures, 
many practical devices and detectors exist in which the 
phonons are expected to move freely only within the 
plane of a membrane, forming a quasi-2D system Q. 
The question how the two-dimensionality of the phonon 
modes influences c-p interaction has been addressed the- 
oretically for certain cases p| [t| [l(| , but no clear exper- 
imental observation of the effect has been reported to 
date, although several attempts have been made pdl[T^|. 

In this paper, we show for the first time experimen- 
tally that the electron-phonon interaction clearly changes 
depending on the dimensionality of the phonons, as ex- 
pected from theory. E-p coupling was measured with the 
help of sensitive NIS tunnel junction thermometry [l3l ]. 
for thin Cu wires on suspended silicon nitride (SiN x ) 
membranes with thickness varying from 30 nm to 750 
nm, which spans the transition from 2D to 3D phonons. 
In addition, samples with identical Cu wires on bulk 
substrates were also measured for comparison. For the 
thinnest membranes, the e-p interaction was strengthened 
in comparison with the bulk samples, and its tempera- 
ture dependence changed significantly, as is predicted by 
the theory @, |, [To|. The change was large enough to 
give indirect evidence that the dispersive (u> ~ fc 2 ), flex- 
ural modes of the membrane likely play a major role in 
the e-p interaction. 

In the presence of stress-free boundaries, the bulk 
transversal and longitudinal phonon modes (with sound 
velocities Ct and q, respectively) couple to each other 
and form a new set of eigenmodes, which in the case 
of a suspended membrane are known as the horizontal 
shear modes (h), an d sy mmetric (s) and antisymmet- 
ric (a) Lamb modes |14| . The frequencies u> for the h 



where fcii is the 



modes are simply ui = Ct + (mir/d) 

wave vector component parallel to the membrane sur- 
faces, d is the membrane thickness and the integer m is 
the branch number. However, the dispersion relations of 
the s and a Lamb modes cannot be given in a closed an- 
alytical form, but have to be calculated numerically. The 
lowest three branches, dominant for thin membranes at 
low temperatures, have low frequency analytical expres- 
sions: u h = c t k\\, lo s = c s k\\, and u a = ^rfcjj, where 

is the effective sound velocity of 
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effective mass for the a-mode "particle". This lowest a- 
mode with its quadratic dispersion is mostly responsible 
for the non-trivial behavior of the e-p interaction 0, [l(| • 
Note that already a single free surface affects the modes 
[HI and the e-p interaction [l6j], as the bulk modes cou- 
ple and form another new set of eigenstates, including 
the surface localized Ray leigh- mode. Thus, the widely 
observed result for e-p power flow P — YjV(T^ — T?) 
from a metal volume V with T e the electron and T p the 
phonon temperature, is not expected to hold even for 
thin enough films on bulk substrates. 

A schematic of the Cu wire samples on suspended sil- 
icon nitride membranes and the used measuring circuit 
is shown in Fig. [1] 17 samples were made on either sus- 
pended membranes or bulk substrates, where nitridized 
(100) Si wafers with 30, 200 and 750 nm thick low-stress 
SiN x top layers were used as the substrate for both cases. 
The suspension of the SiN^ membranes (size 600x300 
/xm ) was achieved by anisotropic backside wet etching of 
the silicon substate in KOH, and the metallic structures 
were fabricated using standard e-beam lithography and 
multi-angle shadow mask evaporation techniques. As the 
e-p interaction strength is sensitive to the thickness and 
disorder level of the metal we minimized its effect 
by evaporating the Cu wires of a specific thickness on 
all the different substrates simultaneously. Ultrathin Cu 
layers (t— 14-30 nm) were used to strengthen the effect of 
the thin membranes. The oxide layer forming the tun- 
nel junction barriers was produced by thermal oxidation 
of Al. Table Q] presents the essential dimensions of the 



2 



samples discussed in this paper, measured by scanning 
electron (SEM) and atomic force (AFM) microscopies. 
The electron mean free path I was determined from the 
resistance of the wire at base temperature 60 mK, using 
the accurately measured dimensions of the wire. 



TABLE I: Parameters for samples. M= suspended SiN x mem- 
brane and B= bulk substrate. B6 had an oxidized Si sub- 
strate. 



Sample 


SiN* d 


Cu t 


V 


I 


r(0.2K) 


r(0.8K) 




(nm) 


(nm) 


KH 3 ] 


(nm) 


(jjs) 


(/is) 


Ml 


30 


14 


2.71 


5.7 


2.6 


0.16 


Bl 


30 


14 


2.46 


4.9 


7.1 


0.030 


M2 


200 


14 


2.44 


4.6 


15.0 


0.11 


B2 


200 


18 


3.67 


4.1 


6.4 


0.045 


M3 


30 


19 


5.50 


11.2 


2.2 


0.30 


B3 


30 


19 


4.62 


9.8 


4.3 


0.034 


M4 


750 


22 


6.09 


10.3 


3.1 


0.030 


B4 


750 


22 


5.87 


8.7 


3.9 


0.013 


M5 


30 


32 


6.09 


22 


1.8 


0.31 


B5 


30 


32 


5.09 


19 


2.7 


0.038 


B6 




32 


7.10 


22 


1.6 


0.031 




FIG. 1: (Color online) A Schematic of the suspended samples 
and the measuring circuit. Red lines are the normal metal 
Cu, light gray Al for SINIS-j unctions and dark gray Al or Nb 
for SN-junctions. 

We used the hot-electron technique [3J] to measure the 
e-p interaction by overheating the electrons by Joule heat 
power P and measuring the resulting electron tempera- 
ture T e . All the samples had two electrically isolated 
Cu normal metal wires next to each other (Fig. 1). The 
longer wire (L — 500/im) was heated by applying a slowly 
ramping voltage across the pair of superconducting Nb 
(or Al) leads in direct metallic contact to Cu, forming 
SN junctions. These junctions provide excellent electri- 
cal, but very poor thermal conductance due to Andreev 
reflection, as the junctions are biased within the super- 
conducting gap A. Thus, due to the lack of outdiffu- 
sion of electrons and the long length of the wire, input 
heat is distributed uniformly in the interior of the wire 
and the electron gas cools dominantly by phonons, in- 
stead of diffusively [l8[ or by thermal photons [l^| . Since 
L » L e _ e , the electron-electron scattering length, elec- 



tron temperature is also well defined without complica- 
tions from non-equilibrium [20j] . In our sample geometry 
the electron temperature is measured with two additional 
Al leads forming a NIS tunnel junctions pair (SLNIS) in 
the middle of heated wire, as a function of input Joule 
power P = IV measured in a four probe configuration. 
The purpose of the short Cu wire, with additional SI- 
NIS thermometer on it, is to give an estimate of the local 
phonon temperature T p , as the e-p power flow depends 
on both T e and T p . 

The current-biased Al SINIS thermometer is ideally 
suited to measure temperature below a few Kelvins, [2| 
due to its high sensitivity (in our DC measurement ~ 0.1 
mK at 0.1 K) and low power dissipation. In addition, 
for all the data here, the SINIS voltage vs. temperature 
response follows the BCS theory without fitting param- 
eters very accurately at least down to ~ 0.2 K, where 
typically saturation sets in. This saturation depends on 
the strength of the e-p interaction (size of thermometer 
and type of substrate) and the amount of filtering, and 
thus we conclude that it is most likely caused by external 
noise heating. For this reason we take the most conser- 
vative approach and assume that all saturation is caused 
by it, in which case we can use BCS theory to convert 
the measured voltage data for all temperatures. 

Even if the electrons lose their energy overwhelmingly 
to the phonons in our sample geometry, it is still pos- 
sible that the measured temperature is not only deter- 
mined by the e-p interaction. This is because the emitted 
phonons could be removed so ineffectively from the mem- 
brane that the phonon transmission becomes a bottleneck 
for the energy flow. Bulk scattering of phonons at low 
temperatures is very weak Q, even for thin disordered 
membranes [Hj], as is boundary resistance for thin films 
on bulk substrates [12, HH. In contrast, almost noth- 
ing quantitative is known about the boundary resistance 
between a thin metal film and a thin 2D membrane, or 
between a thin 2D membrane and a bulk substrate. How- 
ever, it seems clear that if the combined metal film and 
membrane thickness is below the thermal wavelength of 
the phonons, the phonon modes in the two materials are 
strongly coupled, leading to an effectively non-existent 
boundary resistance. Hence, if we check that the mem- 
brane temperature T p is not too high compared to T e 
(effective enough hot phonon removal), we can be confi- 
dent that the measured T e reflects the e-p interaction. 

Figure [2] shows the main result of the measurements, 
with T e and T p plotted vs. the heating power density 
p = P/V for all membrane thicknesses (30 nm, 200 nm 
and 750 nm). In addition, data from a few represen- 
tative bulk samples are shown. Compared to the cor- 
responding bulk substrate sample (B4), T e of the 750 
nm membrane (M4) shows no difference at all, and it 
effectively behaves as bulk. This is reasonable, because 
for the 750 nm membrane the estimated dimensional- 
ity cross-over temperature [H], HH] T cr = Hc t /(2kBd) is 
~ 30 mK, with c t = 6200 m/s for SiN. The phonon 
temperatures T pi however, show a big difference: The 
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FIG. 2: (Color online) Measured electron and phonon tem- 
peratures T e and T p versus the applied heating power density 
in log-log-scale. 
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bulk samples show almost no response from the satura- 
tion value of the thermometer ~ 190 mK, whereas the 
membrane phonons heat up measurably, most likely due 
to the boundary resistance between the membrane and 
the bulk. Nevertheless, this increase in T p for all sam- 
ples is small enough not to influence the e-p interaction. 
For the 200 nm thick membrane (M2) (T cr ~ 110 mK), 
at low heating power densities [p < 40 pW/(pm) 3 ] the 
temperature dependence follows the behavior of the bulk 
sample (B2), although with a difference in the absolute 
value. This shows that the strength of the e-p coupling 
weakens compared to the bulk. At higher powers and 
temperatures (p > 40 pW/(/im) 3 , where T e > 0.6 K), 
T e starts to increase more rapidly in the membrane sam- 
ple, most likely due to the boundary resistance effects. 
The phonons in the 30 nm thick membrane sample (Ml) 
are expected to be in the 2D limit at low temperatures 
(T cr ~ 0.5K), and a clear sign of this can be seen in 
Fig. [5] as a strongly different behavior of the measured 
T e vs. p curve with respect to all other samples. Below 
~ 6 pW/(/im) 3 the e-p coupling is notably stronger (T e 
lower) than in the corresponding bulk (Bl) or any other 
sample, but again at highest temperatures the influence 
of other effects starts to dominate over the e-p coupling. 

To study the temperature dependence of the data in 
Fig. [2] more accurately, we plot the logarithmic deriva- 
tives d(logp)/d(logT e ) in Fig. [3] (a)-(c). For low heat- 
ing powers (T™ >> T") P e - p w T™, where n is the 
power law of the e-p interaction, thus in that regime 
d(\ogp) /d(logT e ) = n. Typically this exponent is n w 5 
for thicker (t > 30 nm) metal films on bulk substrates 
3] , yk [T n| , if the disorder in the film is not too strong 
M rHzUii- From Fig. [3] (a) we first of all see that 
for the 30 nm membrane sample Ml, the difference to 
the bulk sample Bl is very clear. The Ml data has a 



FIG. 3: (Color online) Numerical logarithmic derivatives of 
the measured data in Fig. [2] (a) T e data for Ml and Bl, (b) 
T e data for M2 and B2, (c) T e data for M4 and B4. 



plateau of n ~ 4.5 between p = 0.1 - 6 pW/ (/mi) 3 , while 
for Bl, n continuously decreases from much higher val- 
ues. Note that the strong increase of d(logp)/d(logT e ) 
below p ~ 0.1 pW/(/jm) 3 is caused by the saturation 
of the T e measurement, and not by the e-p interaction. 
The point where n starts deviating from n = 4.5 cor- 
responds to T e « 0.4 K, which is surprisingly consistent 
with the estimated T cr ~ 0.5 K. In contrast, the tempera- 
ture dependence of the 200 nm membrane (M2) and bulk 
(B2) samples [Fig. |3j (b)] are identical with each other 
and with the 30 nm bulk sample (Bl), as long as the 
e-p interaction is dominant (up to 40 pW/(/jm) 3 ). The 
750 nm membrane (M4) and bulk (B4) samples also give 
identical values of n [Fig. [3] (c)]. The difference between 
sample pairs M4,B4 and M2,B2 is caused by the Cu wire 
thickness, which is expected to influence the temperature 
dependence strongly [16l [27j ■ 

Finally, we discuss the effect of the Cu wire thickness 
on the measured e-p interaction. The results for the 
thinnest 30 nm membrane samples, with Cu thickness 
t = 14,19 and 32 nm are shown in Figs [4] (a) and (c). 
It is apparent that the metal film thickness has only a 
minor effect on the e-p interaction on thin membranes, 
and only influences the boundary resistance in the 3D 
limit, by increasing its effect for thicker t, as expected. 
However, for wires on bulk substrates, Figs[4](b) and (d), 
the effect of the Cu wire thickness on e-p interaction is 
more profound. The thinner the Cu film, the more its 
temperature dependence deviates from n = 5, which, for 
comparison, is observed for a more typical t — 32 nm 
Cu wire on oxidized Si (B6). This behavior is qualita- 
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FIG. 4: (Color online) (a) T e versus p — P/V for 30 nm mem- 
brane samples M1,M3,M5. (b) T e versus p for bulk samples, 
from top to bottom Bl (top), B3, B5 and B6 (bottom), (c) 
d(logp)/d(logT) of the data in (a), (d) d(logp)/d(logT) of 
the data in (b). From top to bottom: green line Bl (top), 
magenta B3, blue B5, Red B6 (bottom). In (d) noise has 
been filtered to help the eye. 

tively consistent with the predicted effect of the surface 
phonon modes [16j . but could also depend on the disor- 
der, as the thickening of the film increases the mean free 
path I (Table |T| and pushes the sample closer to the clean 



limit. An apparent exponent as high as ~ 7 could pos- 
sibly be explained by the combination of strong disorder 
and surface modes, but again, detailed theory is lacking. 

In conclusion, we have obtained the first clear evidence 
that the electron-phonon interaction at low tempera- 
tures changes quite significantly when the phonon modes 
become two-dimensional. To quantify the effects, the 
electron thermal relaxation times r = -fVT e / (dP / dT e ) , 
where 7 = 100 J/K 2 m 3 for Cu, are presented in Table U 
for all the samples at two temperatures T e = 0.2 and 0.8 
K. At T e < 0.5 K, the thinnest membranes can have a 
a factor 2-3 strengthening effect, whereas at higher tem- 
peratures the thermal relaxation from membranes can be 
an order of magnitude weaker compared to bulk samples. 
The membrane close to transition region (<i=200 nm) was 
shown to have a weaker (~ factor of two) e-p interaction 
strength than the bulk samples. Thinning the metal film 
on bulk substrates also leads to a sizeable weakening of 
the e-p interaction. The observed power law exponent 
for the 2D limit is consistent with n « 4.5, and is much 
smaller than the corresponding bulk exponent n = 6. .7. 
A reduction by more than a factor one gives indirect evi- 
dence of the importance of the flexural, dispersive Lamb- 
modes for the membrane electron-phonon interaction, in 
agreement with theory [9|, [l0| . 
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